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Abstract 

The Hell— >^HeI recombination of primordial helium plasma {z = 1500 — 3000) is considered 
in terms of the standard cosmological model. This process affects the formation of cos- 
mic microwave background anisotropy and spectral distortions. We investigate the effect 
of neutral hydrogen on the Hell— >HeI recombination kinetics with partial and complete 
redistributions of radiation in frequency in the Hel resonance lines. It is shown that to 
properly compute the Hell— >HeI recombination kinetics, one should take into account not 
only the wings in the absorption and emission profiles of the Hel resonance lines, but also 
the mechanism of the redistribution of resonance photons in frequency. Thus, for exam- 
ple, the relative difference in the numbers of free electrons for the model using Doppler 
absorption and emission profiles and the model using a partial redistribution in frequency 
is 1 - 1.3% for the epoch z = 1770 — 1920. The relative difference in the numbers of free 
electrons for the model using a partial redistribution in frequency and the model using a 
complete redistribution in frequency is 1 - 3.8% for the epoch z = 1750 — 2350. 
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1 Introduction 



The recombination of primordial plasma is a process that ultimately leads to the 
formation of neutral atoms from ions and free electrons due to the decrease in temperature 
through cosmological expansion. This process has three distinct epochs at which the 
fraction of free electrons changes significantly: (1) Helll— >HeII recombination {z ~ 5000 — 
7000), (2) Hell^Hel recombination {z ~ 1500 - 3000), and (3) HII^HI recombination 
(z ~ 900 — 1600), where z is the cosmological redshift. Since other nuclides (D, ^He, Li, B, 
etc.) in the primordial plasma are much fewer in number than and ^He (< 10"^), the 
recombination of hydrogen-helium plasma is usually considered (Zeldovich et al. 1968; 
Peebles 1968; Matsuda et al. 1969). The recombination of other elements is considered 
in isolated cases for special problems, such as the effect of lithium recombination on the 
cosmic microwave background (CMB) anisotropy (Standi et al. (2002) and references 
therein), the formation of primordial molecules (Galli and Palla (2002) and references 
therein), etc. 

The recombination of primordial plasma affects significantly the growth of gravita- 
tional instability and the formation of CMB spectral distortions and anisotropy (Peebles 
1965; Dubrovich 1975). The appearance of the first experimental data on CMB anisotropy 
(Relikt 1, COBE) rekindled interest in the recombination of primordial plasma in the mid- 
1980s. A number of improvements in the model of hydrogen plasma recombination were 
suggested (Jones and Wyse 1985; Grachev and Dubrovich 1991). 

Significant progress in CMB anisotropy observations achieved in the second half 
of the 1990s (BOOMERANG, WMAP) necessitated including a number of subtle effects 
that could affect the recombination of primordial hydrogen and helium at a level of 0.1 
- 1% (Leung et al. 2004; Dubrovich and Grachev 2005; Novosyadlyj 2006; Burgin et al. 
2006; Kholupenko and Ivanchik 2006; Wong and Scott 2007; Chluba and Sunyaev 2006, 
2007, 2008a, 2008b; Hirata and Switzer 2008; Sunyaev and Chluba 2008; Hirata 2008; 
Grachev and Dubrovich 2008). 

One of the most important (for the primordial plasma recombination kinetics) ef- 
fects considered in recent years is the absorption of Hel resonance photons by neutral 
hydrogen, which leads to an acceleration of the Hell— s^Hel recombination (Kholupenko 
et al. 2007 [hereinafter KhIV07]; Switzer and Hirata 2008; Rubino-Martin et al. 2007). 
Recently, this effect was taken into account by Wong et al. (2008) in the recfast compu- 
tational code developed by Seager et al. (1999). This code is most widely used to compute 
the primordial plasma recombination kinetics when the CMB anisotropy is analyzed. To 
take into account the effect of neutral hydrogen on the Hell— >HeI recombination kinetics, 
Wong et al. (2008) used a simple approximation formula with adjustable parameters. 
The recfast modified in this way allows the results of computations with the multilevel 
CO d^ by Switzer and Hirata (2008) to be quickly and accurately reproduced for any rea- 
sonable values of the cosmological parameters. Nevertheless, the approach by Wong et 
al. (2008) is inapplicable for more detailed studies of the helium recombination, since 
their formula is not universal for all resonance transitions in Hel, but can be used only 
in describing the absorption of Hel 2^P — > l^S resonance photons by neutral hydrogen. 
When using the multilevel codes by Switzer and Hirata (2008) and Rubino-Martin et 
al. (2007), which allow the absorption of Hel nP ^ IS (n > 2) resonance photons by 



^By the multilevel code we mean a computational program that uses a multilevelmodel atom. 
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neutral hydrogen to be taken into account, much of the radiative transfer calculations in 
the Hel resonance lines are performed numerically, which is computationally demanding 
and time-consuming. These circumstances forced us to seek for an analytic approach to 
the problem of including the effect of neutral hydrogen on the Hell— >^HeI recombination 
kinetics that, on the one hand, would be more universal than the approach of Wong et al. 
(2008) (i.e., would allow the effect of the absorption of Hel nP IS resonance photons 
(where n > 2, and not only n = 2) on the Hell^Hel recombination to be estimated) and, 
on the other hand, would not reduce the speed and accuracy of the primordial plasma 
recombination computations for various parameters of the cosmological model. This ap- 
proach was implemented on the basis of the papers by Chugai (1987) and Grachev (1988), 
who analytically investigated the diffusion of resonance radiation in the presence of con- 
tinuum absorption. In this paper, we present extended and more detailed justifications 
of the key suggestions made in KhIV07. We take into account the fact that the scatter- 
ing in the Hel resonance lines (for transitions in the singlet structure of the Hel atom) 
occurs with a partial redistribution in frequency, which turns out to be important for the 
results (Switzer and Hirata 2008; Rubino-Martin et al. 2007). Thus, our goal is to nu- 
merically compute the Hell— s^Hel recombination kinetics using analytic formulas (Chugai 
1987; Grachev 1988; this paper) to allow for the peculiarities of the radiative transfer in 
the Hel resonance lines (a partial redistribution of Hel resonance photons in frequency 
and their absorption in the neutral hydrogen continuum). 

2 Physical model of Hell^Hel recombination 

In the process of cosmological recombination, the plasma deviates from its ionization 
equilibrium. A necessary condition for this deviation is plasma opacity for the intrinsic 
resonance recombination radiation. This means that the emitted resonance photon is 
absorbed by another neutral atom almost instantly (compared to the characteristic re- 
combination and ionization time scales). Since the recombination radiation is excessive 
with respect to the equilibrium background with a blackbody spectrum, the populations 
of excited atomic states exceed their equilibrium values. The excess of the excited-state 
populations entails an increase of ionization fraction compared to its equilibrium value 
and, accordingly, leads to a delay of recombination. In this situation, the plasma recom- 
bination can no longer be described by the Saha formula and kinetic equations should be 
invoked to describe the behavior of the excited-state populations of neutral atoms and the 
plasma ionization fraction. The total recombination rate Jtot [cm^'^s^^] (dependent on z) 
is determined by the sum of the recombination rates (given the forward and backward 
reactions) to all of the bound Hel atomic states. 



The atomic transition rate to the ground state J_^i is defined by the sum of the 
transition rates from all of the excited states and the continuum. 
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Since the number of Hel atoms in excited states at the Hell— >HeI recombination epoch is 
not accumulated (no more than 10~^ of the Hel atoms are in each excited state), it may 
be concluded that Jtot — J^i- The quantities J„i in sum (2) are the differences between 
the direct, J„^i, and reverse, Jn^i, transition rates, 

Jnl = Jri^l — Jn^l (3) 

The estimates made by Zeldovicli et al. (1968) and Peebles (1968) and subse- 
quently confirmed by numerical calculations using multilevel model atoms (Grachev and 
Dubrovich 1991; Seager et al. 2000) showed that a simplified recombination model (the 
so-called three-level model; Zeldovich et al. 1968; Peebles 1968; Matsuda et al. 1969; 
Seager et al. 1999) could be used to calculate the ionization fraction as a function of 
time. In this model, the recombination rate for helium (the model energy level diagram is 
presented in Fig. 1) is determined by the following processes: the two-photon 2^5' — > 1^5* 
transitions and the one-photon 2^P — > V-S and 2^P — > 1^5", i.e., three terms remain in 
sum (2) and the following formula is valid: 

Jtot — Jag + Jbg + Jh'g (4) 

where the subscripts denote the following states: g = l^S, a = 2^S, b = 2^P, b' = 2^P 
(see Fig. 1). 

A proper allowance for the resonance transitions requires a joint analysis of the 
kinetics of transitions and radiative transfer in the 2^P 1^5* and 2"^P 1^5* lines by 
including a number of peculiar factors, with the cosmological expansion and the absorption 
of Hel resonance photons by neutral hydrogen (HI) being the most important of them. 

According to Eq. (3), the two-photon Hel 2^5" — > l^^S" transition rate can be 
calculated using the formula (for convenience, the common factor was taken out of the 
brackets): 

Jag = Aa, (^Na - fVagNg^ (5) 

where Aag = 51.3 c^^ is the coefficient of the spontaneous two-photon 2^S l^S decays, 
Na [cm~^] is the population of the 2^S state, Ng [cm^^] is the population of the l^^S" state, 
Qa — 1 is the statistical weight of the 2^5" state, g^^, = 1 is the statistical weight of the 
1^5" state, rjag is the equilibrium photon occupation number at the 2^5" l^S {a — > g) 
transition frequency. 

For optically thick transitions (for Hel at the Hell^Hel recombination epoch, these 
are the n^P — > l^S and 2^P — > 1^5" transitions and, when the three-level model is used, 
the b ^ g and b' ^ g transitions, respectively) the rates of the direct and reverse processes 
appearing in Eq. (3) are very close (their relative difference can reach 10~^, depending 
on n and the instant of time under consideration), because the occupation numbers of 
the photon field in Hel lines (including both equilibrium and intrinsic Hel recombination 
radiations) are close to their quasi-equilibrium values, rjfg ~ N^/ gnNi. 

Thus, if J„i are calculated from Eq. (3), then two relatively close numbers often 
has to be subtracted. In this case, a significant loss of the computation accuracy is possible 
(Burgin 2003). Therefore, the following formula that is devoid of the above shortcoming 
is used to consider the kinetics of such optically thick transitions: 

Jf, = PfAf, (^/ - 1^/.^^) (6) 
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where the subscript / denotes the nP state. Afg is the Einstein coefficient for the spon- 
taneous f ^ g transition, Nf is the population of state /, and gj is the statistical weight 
of state /. The quantity Pfg is the probability of the uncompensated f ^ g transitions. 

The quantity Pfg can be calculated by jointly considering the radiative transfer 
equation in the f ^ g line and the balance equation for levels / and g. It contains 
information about the effect of the intrinsic resonance plasma radiation on the transition 
kinetics with allowance made for the circumstances that accompany the radiative transfer, 
such as the cosmological expansion, the absorption of Hel resonance photons by neutral 
hydrogen, etc. 

Since the / — > 5' transition rate Jfg, which, in turn, determines the recombination 
rate Jtot, directly depends on Pfg, a proper calculation of Pfg for the primordial plasma 
conditions is one of the most important subgoals of the cosmological recombination theory. 

The methods for including Pfg in the multilevel code that computes the kinetic 
equations for the full system of levels can be found in Seager et al. (2000), Switzer and 
Hirata (2008), and Rubino-Martin et al. (2007). The kinetic equation that describes the 
Hell— >HeI recombination in terms of the simplified model and the methods for including 
Pfg in it can be found in KhIV07 and Wong et al. (2008). 



3 Kinetics of Hel 2P^1S resonance transitions 

Let us consider the Hel f ^ g resonance transition kinetics using Hel 2P^1S (i.e. b ^ g 
and b' ^ g), which mainly determine the recombination rate Jtot, as an example. A joint 
analysis of the balance equation for levels / and g and the radiative transfer equation for 
the Hel f ^ g line leads to the following formula for the probability of the uncompensated 
f ^ g transitions: 

Pf9= / —irz — -^^f9{^)d^ (7) 

Jo K'H + K'HeJ 

where kh [cm~^s~^Hz^^] is the absorption coefficient of photons {hu > 13.6 eV) by neutral 
hydrogen during ionization, KhcJ is the absorption coefficient of photons by neutral helium 
in the f ^ g line, r/ is the optical depth for the absorption of Hel resonance photons 
(including the absorption by both helium and hydrogen atoms), and '0/ff(^) is the emission 
profile in the Hel / — > (? line (/ ipfg{u)di' = 1). 

The coefficient is given by the formula 

87rj/^ /NAT / \ 

K'H = —^(yH{y)NH,lS (8) 

where cr//(z/) is the photoionization cross section of the HI ground state by a photon with 
frequency v, Nh,is is the number density of hydrogen atoms in the ground state, which 
is equal, with a good accuracy, to the neutral hydrogen number density Nhi- 
The coefficient KhcJ {f = b or b') is given by the formula: 

f^HeJ = ^AfgNg<l>fg{iy) (9) 

99 

where (pfg{i^) is the absorption profile in the Hel f ^ g line (/ (j)fg{i')dh' = 1). The 
ground-state population Ng is equal, with a good accuracy, to the total neutral helium 
number density NhcI- 
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The optical depth Tf is given by the formula 

where H{z) = HQ^JVtx + + zf + fi^e/(l + z)"^ is the Hubble constant (the parame- 
ters Hq, ^a, ^rei are described below in Table Ej), and the parameter z' is defined 
by the equality z' = {1 + z)v' jv — 1. 

For the convenience of the subsequent consideration, let us represent Pfg as the 
sum of two terms: Pfg = + Pfg'^, where Pj^ and PJ^'^ are given by the formulas: 







The quantity P^ is the mean (averaged over the profile ^Jfgiy)) destruction proba- 
bility of a Hel f ^ g resonance photon as a result of its interaction with neutral hydrogen. 

The quantity PJg'^ is the escape probability of Hel f g photons from the line 
profile due to the cosmological expansion. Note that the definition of this quantity differs 
from the classical definition of the Sobolev photon escape probability from the line profile 
Pfg (Rybicki and dell Antonio 1993; Seager et al. 2000), being its generalization to 
the case that includes the absorption of Hel resonance photons by neutral hydrogerlf] 
(in this sense, PJg'^ may be called a modified photon escape probability from the line 
profile). If the absorption of Hel resonance photons by neutral hydrogen is negligible (i.e., 
f^ul f^Hej ^ 1), then the formula for PJ^'^ takes the classical form: 

P;fU,=o= / e^v{-rHeAi^))i^fgWv^P^g^THlf{l-e^v{-rHej)) (13) 

where Tneji^) is the optical depth for photon absorption in the Hel f ^ g line as a 
function of the frequency and ThsJ is the total optical depth for photon absorption in the 
Hel f ^ g line given by the formula: 

Note that Switzer and Hirata (2008) and Rubino-Martin et al. (2007) used a 
different splitting into terms, namely, Pfg = Pfg + APj^^, where ^Pfg'^ is the correction 
to the Sobolev escape probability due to the presence of continuum absorption. 

3.1 The Modified Photon Escape Probabihty from the Line Pro- 
file 

Using ([9]) and (fTOl l. we can transform Eq. ( fT2i ) to 

Pjf = r ( _J^JtJ_) ' ( tlM) exp (-r') dr' (16) 



^Note that Pjg'^ is proportional to (1 — I^), where II is the integral introduced by Switzer and Hirata 
(2008). 



This integral can be roughly estimated from the formula 



Pff = (1 + I'T'^hIj (1 - i-THej)) 



(16) 



where the parameter 7 is the ratio of the helium and hydrogen absorption coefficients at 
the central frequency of the f ^ g line and is given by the expression: 



The z dependence of 7 for the 2^P l^S and 2^P — > l^S transitions is presented in Fig. 



If the amount of neutral hydrogen is negligible (so that 7 -C 1) then Eq. (fT6l) . as 
has been noted above, transforms into the standard expression for the Sobolev photon 
escape probability from the line profile - ffT3l) . If there is much neutral hydrogen 
(so that 7 > 1), then the value of PJg'^ is lower than the typical Sobolev value of P^, 

calculated from In KhIVOT, the factor (1 + 7"^) in Eq. ^ was discarded 

(i.e., in fact, Eq. (fT3|l )was used). For the Hel 2^P l^S transition this neglect is valid 
throughout the Hell— >HeI recombination epoch, because the values of 7 for this transition 
are large (see Fig. [2]). For the Hel 2^P l^S transition, this neglect is valid at the early 
Hell— >^HeI recombination epoch (z > 1900), when the values of 7 for this transition are 
large (> 10, see Fig. [2|). At the late Hell— >HeI recombination epoch {z < 1750), when 
the values of 7 for the Hel 2^P l^S transition are small (< 1, see Fig. [2|), the classical 
expression for the probability of the uncompensated transitions (fT3l) is inapplicable for 
the Hel 2^P l^S transitions. At the same time, however, the 2^P l^S transition 
rate Jbg, which directly depends on the product AbgPj^, is so large that it completely 
determines the recombination rate Jtot- Therefore, the accuracy of calculating P^,'g'^ (and, 
accordingly, the 2^P l^S transition rate Jb'g) does not play significant role. 

3.2 The Destruction Probability of a Hel Resonance Photon dur- 
ing Its Interaction with Neutral Hydrogen 

The integral expression (fTTl) for P^ can be transformed to the following form convenient 
for both numerical and approximate analytical integrations: 



A further refinement of the form of the function P^ i^) depends on the approximation 
in which the absorption profile (f)fg{i^) is taken into account and, even more importantly, 
on the specific form of the emission profile ipfg{i^) determined by the physical conditions 
under which the radiation is scattered and transferred in the f ^ g resonance line. 

3.2.1 Complete redistribution in frequency: The Doppler profile 

In the case of scattering with a complete redistribution in frequency, if the Doppler profile 
is used as the absorption, 0/g(z/), and emission, i^fgii'), profiles (as was done in KhIVOT)), 
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7 = 



knei^fg) ^ (gf/gg) AfgNHel(l)fg{jyfg)(^ 

knii^fg) aH{vfg)Mv)gNHi 



(17) 



El 




(18) 



the expression for Pf takes the form: 



/'o"(7)W"" ;;"^^T^^ y (19) 

i_oo l + 7exp(-y2) 

where the subscript D stands for "Doppler", while the subscripts / and g were omitted, 
since the function (7) is universal for all resonance transitions when the Doppler profile 
is used. This function can be approximated by the expression 

pE{i) = {i+Pi'r' (20) 

where the parameters p, q depend on the 7 range. Their values are given in Tabled! 



Table 1: Parameters of the approximation of Pj 



Range 


of 7 




P 


q 


< 7 


< 5 - 102 




0.66 


0.9 


5-102 


< 7 < 5 • 


■ 10^ 


0.515 


0.94 


5-10^ 


< 7 < 5 • 


■ 10^ 


0.416 


0.96 


5-10^ 


< 7 




0.36 


0.97 



The asymptotics of P§ (7) for 7 ^ 00 is given by the expression (see, e.g., Ivanov 

1969): 

Pf(7)^^7-'Vl^ (21) 

The results of the calculations of the function P^ij) for various 7 ranges are 
presented in Figs. [3] and [H 



3.2.2 Complete redistribution in frequency: The Voigt profile 

In the case of scattering with a complete redistribution in frequency, if the Voigt profile 
is used as the absorption, (pfgiy) and emission, ipfgiy) profiles, the probability P^ can be 
approximately calculated from the formula 

P^^Pg + P^ (22) 

where the subscript V stands for "Voigt". The quantity P^ is attributable to the inclusion 
of the Voigt profile wings (the subscript W stands for "wings") and is given by the formula 




where a is the Voigt parameter, which is defined by the ratio of the natural (the quantum- 
mechanical damping constant defining the mean level lifetime) and Doppler widths of level 
f: a = Tf/inAuDj, ai ~ 1.6 is an adjustable parameter whose value is chosen from the 
condition for the best agreement between the values of P^ obtained from Eqs. (fTSll (by 
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numerical integration) and (l22l) . Equation ([23l) can be derived by taking into account 
the fact that almost the entire contribution to 4>fg{^) ^it | (z^ — i^fg) / Au^jl > | ln(a)| is 
provided by the Lorentz wings, i.e., the following approximation is valid: 

y) ~ ^ (^y _ y^g) (24) 



TT 



The thermal width Au^j in Eq. (l24l) is given by the expression (Lang 1978): 



A.z„ = ^J^ + |k^ (25) 

C V "^He 3 

where mue is the helium atomic mass, T is the temperature of the medium, and Vt 
is the root-meansquare turbulent velocity (if the distribution of turbulent velocities is 
Maxwellian). In our calculations, we assumed that = 0. 

The asymptotics of Py (7) for 7 00 is given by the expression (see, e.g., Ivanov 

1969): 

~ T^^'^a'l^r"^ (26) 

The results of the calculations of the function Py (7) for the Hel 2^P — > l^S and 2^P — > 
1^5* and transitions (and various 7 ranges) are presented in Figs. [3] and [4] respectively. 

3.2.3 Partial redistribution in frequency: The Wong-Moss-Scott approxima- 
tion 

Switzer and Hirata (2008) and Rubino-Martin et al. (2007) showed that a partial (rather 
than complete) redistribution in frequency occurred at the Hell— >HeI recombination epoch 
in the Hel n^P lines. These authors performed a significant fraction of their calcu- 

lations (in particular, the radiative transfer calculations) numerically. Since the problem 
under consideration is complex, this is computationally demanding and time-consuming 
(it takes about one day to compute the radiative transfer in Hel lines for one cosmological 
model; Rubino-Martin et al. 2007). This approach is too resource-intensive to be used in 
the three-level recombination model incorporated, in particular, in recfast (Seager et al. 
1999; Wong et al. 2008), from which a high speed of computations (no more than a few 
minutes per cosmological model) is demanded at the required accuracy of about 0.1%. 
This necessitates seeking for an analytic solution to the problem of radiative transfer in 
a resonance line in an expanding medium in the presence of continuum absorption or at 
least a suitable approximation that would describe satisfactorily the results of computa- 
tions with multilevel codes (Switzer and Hirata 2008; Rubino-Martin et al. 2007). This 
approximation was found by Wong et al. (2008): based on an approximation formula of 
form (I20D, Wong et al. (2008) found that at p = 0.36 and q = 0.86 for the probability P^ 
of the uncompensated 2^P l^S transitions (let us denote it by P^ms^ where the sub- 
script WMS stands for ' Wong-Moss-Scott"), the simplified (three-level) model describes 
well the results of computations with the multilevel code by Switzer and Hirata (2008) for 
any reasonable cosmological parameters. The values of p = 0.66 and q = 0.9 (KhIV07) 
were used to determine the 2^P ^ l^S transition probability. 

The asymptotics of Pwms^i) for 7 — 00 is given by the expression: 

Kms = 0.367-°-«^ (27) 
The results of the calculations of the function Pwms ^^^^ presented in Fig. [3l 
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3.2.4 Partial redistribution in frequency in the wings: The Chugai-Grachev 
approximation 

The Wong-Moss-Scott approximation cannot be used to investigate the kinetics of Hel 
n^P — > 1^5* resonance transitions for n > 3, since the parameters p and q in Eq. (I20l ) 
should have different values unique for each specific n in this case. This necessitates 
seeking for an expression for the probability of the uncompensated transitions based on 
an analytic solution of the radiative transfer equation in the Hel resonance lines at the 
Hell— >^HeI recombination epoch. 

Being formulated in full, this problem is very complex, since it requires including 
a large number of processes, such as the Hubble expansion, the partial redistribution in 
frequency due to the thermal motion of atoms, the Raman scattering, the continuum 
absorption, the recoil upon scattering, etc. In this case, the problem requires considering 
an integrodifferential equation in which the redistribution function cannot be expressed in 
terms of elementary functions. Therefore, in this formulation, it is solved mainly through 
computer simulations (Switzer and Hirata 2008; Rubino-Martin et al. 2007). Neverthe- 
less, a number of simplifications and approximations make it possible to reformulate the 
problem in a form that allows a completely analytic solution. Thus, for example, Chugai 
(1987) and Grachev (1988) considered the diffusion of resonance line radiation in the 
presence of continuum absorption. In comparison with the complete formulation of the 
problem of radiative transfer in a line, these authors disregarded the following effects: (1) 
the Raman scattering was disregarded; (2) a differential expression derived in the diffusion 
approximation was used instead of the exact integral term describing the redistribution 
of photons in frequency due to the thermal motion of atoms (see, e.g., Varshalovich and 
Sunyaev 1968; Nagirner 2001), i.e., the final equation has the form of a frequency dif- 
fusion equation for photons (Harrington 1973; Basko 1978); (3) the Doppler core of the 
Voigt absorption profile was fitted by a delta function; and (4) when formulating the 
mathematical model, Grachev (1988) disregarded the expansion of the medium, although 
Chugai (1987) previously took it into account. Comparison of the papers by Chugai 
(1987) and Grachev (1988) in this aspect shows that the contributions to the probability 
of the uncompensated transitions from the expansion of the medium and the continuum 
absorption can be approximately taken into account independently of each other. Chugai 
(1987) derived the following formula for the probability of the uncompensated transitions 
(note that Chugai (1987) and Grachev (1988) used notations differing from each other 
and from the notation of this paper): 

P^ = 1.217ai/S"'^^ (28) 

where the subscript C stands for "Chugai". The functional dependence in Eq. (l28ll was de- 
termined from qualitative considerations, while the numerical coefficient was determined 
by numerically solving the diffusion equation. 

Grachev (1988) analytically derived the formula: 

= (8Aa)VV-/«7-^/^i^ 2, . + 2, (29) 
where the subscript G stands for "Grachev", A is the single-scattering albedo, F is the 
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hypergeometric function, T is the gamma functioE@, and the parameter s is defined by 
the formula: 




2) at large 7 (7 > 71) 



-1 



(34) 



We see that at small 7 (7 ^ 71, which corresponds to 7 ^ 5 ■ 10^ for the Hel 
2^P l^S transition) and A close to unity, Eq. f l33l ) matches the solution obtained by 
Chugai (1987). 

It should be noted that Eqs. (1281) and f l29l ) were derived by Chugai (1987) and 
Grachev (1988) for 7 at which 07 1 and cannot be applied at 07 < 1. This can be 
seen from the following: 1) when 7^0, and Pq tend to infinity; 2) when the Voigt 
parameter tends to zero (a — 0), Pq and Pq also tend to zero, while the probability P^ 
that they describe tends to that corresponds to the correct description of the Doppler 
core (in terms of the Doppler profile rather than the delta function, as was done by Chugai 
(1987) and Grachev (1988)). It should also be noted that the correct description of the 
Doppler core leads to a correction of about 10% to Pq for the Hel 2^P — > l^S transition 
in the range 10^ < 7 < 10^. 

The results of the calculations of the function Pq , its asymptotics (1331) and f l34l) . 
and the relative error in Pq calculated from the approximation formulas f l30l) and (I3T1 ) 
compared to the calculations based on the exact formula (l29l) are presented in Fig. [5l 
Comparison of asymptotics (l33l) and (l34l) (Fig. [5l) shows that the 7 dependence of Pq , 
changes significantly at 7 close to 71 (e.g., on a logarithmic scale, the slope d hi Pq / d hi 'j 
changes from -0.75 to -1). This should be taken into account when the Chugai-Grachev 
approximation is used to calculate the Hell— >^HeI recombination kinetics, because 7 (for 
the Hel 2^P l^S transition) takes on values close to 71 ~ 5 • 10* during the Hell— >HeI 
recombination (at the epochs z = 1900 — 2100) (see Fig. [2|). 

3.2.5 Partial redistribution in frequency: An approximate allowance for the 
Raman scattering 

Since Chugai (1987) and Grachev (1988) (1) disregarded the Raman scattering and (2) 
described the central region of the absorption profile by ^-function, Eq. (l29ll gives an 

''The following approximation formulas expressible only in terms of elementary functions can be used 
to calculate the expressions containing special functions in Eq. ((29l) : 



These formulas are valid in the interval of s [—1/4; 00) with a relative accuracy of at least 0.5%. 



r(s + 3/2) 

r(s + 2) 



(s + 1.28) 



1/2 



(30) 




(31) 
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underestimated value compared to the probability , determined when solving the full 
problem, in which the Raman scattering contributing to the formation of broader wings in 
the emission profile ipfgi^u), is taken into account and the central region of the absorption 
profile (Doppler core) is described by the Voigt profile. To estimate the contribution 
from the Raman scattering and the Doppler core to the probability of the uncompensated 
transitions P^ , let us represent Eq. (fTSll as: 

P^c^PE + PS + Pi (35) 

where P^, Pq and P^ are the contributions to the probability P^, from various physical 
effects. These quantities are explained in detail below. 

1) The contribution P^ can be defined by the formula: 

PE= I , ^ I. f w^I^ I ^^ dv , A ~ [ufg - 3Aud; Ufg + 3Aud] (36) 

7^1+(0/g(Z/)/0^g(Z/^3))7 

where A is the Doppler core region. The quantity P^ includes the contribution from the 
absorption of Hel resonance photons from the Doppler core region to the probability of the 
uncompensated Hel 2^P — > 1^5, which is disregarded in the Chugai-Grachev approach. 
The quantity P^ can be calculated using Eq. ( l20l l. 

2) The contribution Pq can be defined by the formula: 

PB = f , , , , f^w^I^ , T^du , B ~ [ufg-ucG-, iyfg-3AiyD]U[ufg+3AuD; ^fg+^CG] 

(37) 

The quantity P^ includes the contribution from the region of the "near" wings P, where 
the redistribution in frequency due to the thermal motion of atoms plays a crucial role 
in forming the emission profile, just as in the Doppler core region. Pq can be calculated 
using Eq. ( l29l l. The characteristic frequency ucg specifies the boundaries of the frequency 
region B in such a way that the redistribution in frequency due to the thermal motion 
of atoms considered by Chugai (1987) and Grachev (1988) has a decisive effect on the 
formation of the emission profile in the frequency range [vjg — vcg] ^fg + ^cg], while the 
Raman scattering has a decisive effect on the formation of the emission profile in the 
frequency ranges [0; Ufg — vcg] and [vfg + vcg] oo]. We calculate vcg below. 

3) The contribution P^ can be defined by the formula: 

'"^ Lr^^wMh^,'" ' C.|0;.,,-.cclU[., + .cc;ool (38) 

The quantity P^ includes the contribution from the region of the "far" wings C, where 
the Raman scattering (the subscript R stands for "Raman") plays a crucial role in forming 
the emission profile. 

Since the emission and true absorptioE@ profiles in the ranges of integration in (l38ll 
are defined by the expressions ipfgiy) = (1 — A)0y(z/), and (pfgiy) = (1 — A)0y(z/) (where 
0y(i/) is the Voigt profile), we obtain 



PR = 2(1 - A)^/^^3^:^ - arctan(/(a, A,7))) (39) 



^The term "true absorption" is used here in the same sense as that in Ivanov (1969) and Nagirner 
(2001). 
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where /(a, A, 7) is defined by the expression 

(40) 



/(a, A, 7) = AiTrV^ ((1 - A)a7)-'/' (^"j 



Here, is the correction factor that takes into account the fact that the two effects 
(both the redistribution in frequency due to the thermal motion of atoms and the Raman 
scattering) give comparable contributions to the probability at frequencies close to the 
boundary frequency Vjg — Vcg (and, accordingly, v^g + Vcg) (i-e., there is no well-defined 
boundary between the zones of infiuence of these effects). 

Using the results by Chugai (1987) and Grachev (1988), one can show that the 
integrand in (l35l) in the central frequency region \vjg — vcg\ ^jg + ^cg\ depends on the 
frequency as 

-expf-( „„.__.,L^,.,j 1 (41) 



where x = [y — Ufg)/ /S.U£i. 

In the frequency ranges [0; Vfg — i'cG] and [I'fg + i'cG'-, 00], where the effect of Raman 
scattering prevails, the integrand in (l35l ) depends on the frequency as 

^fai^) a{l-\) 

l + {(pfg{T^)/<Pfg{vfg))^^ ^x^ + a{l-\)-i^ ^ ' 

Comparing Eqs. (|4T1) and (l42l) . we can find the characteristic frequency ucg in the form 

ucG - 2'/\-^/\\a-iY'^^/]^^UD (43) 
Substituting (l43l) into (l40ll . we obtain the final expression for /(a, A, 7): 

/(a, A, 7) ~ Ai2i/Vi/8(1 - A)-i/2ai/4 (a7)-i/^ (44) 

The coefficient Ai = 0.5 can be determined from the condition for the best agree- 
ment between the dependences Pbg{z) = {P^g'^{z) + P^{z)) (see Fig. [7]), calculated here 
and in Rubino-Martin et al. (2007). 

It should be noted that we assume here that if the fraction of coherent scatterings 
(i.e., the singlescattering albedo) is equal to A, then the fraction of Raman scatterings 
is equal to (1 — A), i.e., only the radiative transitions are taken into account, while the 
effect of the transitions produced by electron collisions is considered negligible. This 
approximation is valid at the Hell— >^HeI recombination epoch. In the general case where 
the effect of electron collisions is not negligible, the fraction of Raman scatterings is not 
equal to (1 — A). This should be taken into account when the above formulas are used. 

As the fraction of coherent scatterings A changes, Eq. (l35l ) has the following 
asymptotics: 

1) As the fraction of coherent scatterings A tends to unity, Eq. (1351) turns into the sum 
of and P^ given by Eq. (l28l) (Chugai 1987). This limit corresponds to the absence 
of Raman scatterings, i.e., the resonance photons are redistributed in frequency solely 
through the thermal motion of atoms. 

2) As the fraction of coherent scatterings A tends to zero, Eq. (I35l l tends to Py given by 
Eq. f l22l ). This limit corresponds to a completely incoherent scattering, which leads to a 
complete redistribution of resonance photons in frequency. As a result, the Voigt emission 
profile is formed. 
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Table 2: Parameters of the standard cosmological model 



Description 


Desi^ 


^nation 


Value 


iotal matter density 


^kot 




1 


(in units of critical density) 








Baryonic matter density 






0.02-0.06 


Nonrelativistic matter density 




= ^CDM + 


0.27 


Relativistic matter density 


^rel ' 


— ^^'Y 1 ^^1/ 


8.23 ■ 10-^ 


Vacuum-like matter density 






0.73 


Hubble constant 


Ho 




70 km/s/Mpc 


CMBR temperature today 


To 




2.726 K 


Helium mass fraction 


Y 




0.24 



4 The Cosmological Model 

We performed all calculations within the framework of standard cosmological ACDM 
models. The calculation results presented in Figs. [6l [7] and [8] were obtained using the 
cosmological parameters from Rubino-Martin et al. (2007) (since these figures reflect, in 
particular, the comparison of our results with those of Rubino-Martin et al. (2007)). The 
results of calculations presented in Fig. [9] were obtained using the cosmological parameters 
from Table H 

5 Results and Discussion 

The calculated destruction probabilities of resonance photons as they interact with neu- 
tral hydrogen are plotted against the ratio of the helium and hydrogen absorption 
coefficients 7 in Figs. [3] and [4] for various models of photon redistribution in frequency. As 
we see from these figures, the probabilities are lowest and highest when the Doppler 
and Voigt profiles, respectively, are used to describe the absorption and emission coef- 
ficients. This is because the Doppler profile "has no wings", while the Voigt profile has 
slowly descending Lorentzian wings. Therefore, the calculated P^ for any other models of 
redistribution (for a partial redistribution) in frequency (including the actual dependence 
P^(7)) should lie between the curves corresponding to the use of the Doppler and Voigt 
profiles (i.e., the inequality P^ < P^ < Py should hold). 

As was shown by Switzer and Hirata (2008) and Rubino-Martin et al. (2007), the 
partial redistribution in frequency should be taken into account when the probability P^ 
is calculated for the Hel 2^P ^ l^S transition, because the fraction of coherent scatterings 
in the total number of resonance photon scatterings in the line is high, A ~ (1 — 2.5-10"^). 
The probability of the uncompensated Hel 2^P l^S transitions can be calculated by 
taking into account the partial redistribution in frequency using Eq. (l35ll or in the Wong- 
Moss-Scott approximation (l27l) (Fig. [3]). In the range of values 10^ < 7 < 10^ important 
for the Hell— i>HeI recombination kinetics, the calculated values of P^ and P^ms ^'^^ 
satisfactory agreement (differ by no more than 50%, although the values of the quantities 
themselves change by two orders of magnitude). 



14 



The scattering in the Hel 2^P ^ l^S resonance line occurs at an almost complete 
redistribution of photons in frequency due to the high fraction of incoherent (Raman) 
scatterings, (1 — A) ~ 1 (and, accordingly, the low fraction of coherent scatterings, A ~ 0). 
Therefore, Eq. (I22l ) which was derived by assuming a complete redistribution of resonance 
photons in frequency, is used to calculate the probability (equal to Py ) for this 
transition. Figure [4] presents the calculated values oi Py and P§ for the Hel 2^P — > l^S 
transition. We see from Fig. [4] that the curves for P^ and Py coincide in the range of 
values 10~^ < 7 < 10"^ important for the transition under consideration at the Hell— i>HeI 
recombination epoch (see Fig. [2|). This is because the contribution from the wings P^r 
to the probability Py at these values of 7 is negligible compared to the contribution 
from the Doppler core P§ . In turn, this is because the Voigt parameter is small for 
the Hel 2^P — > l^S* transition: a ~ 10~^. The contribution from the wings P^ to the 
probability P^ begins to have an affect at such values of 7 that 07 > 1 (this is true not 
only for a complete redistribution in frequency, but also for a partial redistribution in 
frequency). Thus, in Hell— >^HeI recombination calculations, the probability P^ for the 
Hel 2^P ^ l^S transition can be calculated with a sufficient accuracy from Eq. (I20l ) (as 
was done in KhIV07), although it would be more appropriate to calculate this probability 
from Eq. (l22ll . 

The main result of this paper is the dependence of the relative number of free elec- 
trons NJN]E on redshift z for the Hell— >^HeI recombination epoch (Fig. [6|), calculated by 
including the effect of neutral hydrogen when using various models for the redistribution 
of Hel resonance photons in frequency upon their scattering in the Hel 2^P l^S line. 
Figure [6] leads us to the following conclusions: 1) The effect of neutral hydrogen on the 
Hell— >^HeI recombination turns out to be significant, since including it changes N^/Nh 
by 1 — 2.1% for the epochs z = 1650 — 1950 compared to the recombination scenario 
in which it is disregarded (Dubrovich and Grachev 2005; Wong and Scott 2007). This 
change in N^/Nh is significant for a proper analysis of the experimental data on CMB 
anisotropy that will be obtained from the Planck experiment scheduled for 2009; 2) The 
calculated z dependence of N^/Nh at the Hell— s>HeI recombination epoch turns out to 
be sensitive to the model for redistribution of Hel resonance photons in frequency. Thus, 
for example, the relative difference in N^/Nh for the model using the Doppler absorption 
and emission profiles (KhIV07) and the model using a partial redistribution in frequency 
(Switzer and Hirata 2008; Rubino-Martin et al. 2007; Wong et al. 2008; this paper) is 
1 — 1.3% for the epochs z = 1770 — 1920. The relative difference in N^^/Nh for the model 
using a partial redistribution in frequency and the model using a complete redistribution 
in frequency for the Hel 2^P l^S resonance transition (Rubino-Martin et al. 2007; 
this paper) is 1 — 3.8% for the epochs z = 1750 — 2350. These differences are significant 
for the analysis of the experimental data from CMB anisotropy measurements during fu- 
ture experiments (Planck etc.). This suggests that simple models for the redistribution 
of resonance photons in frequency (such as those using the Doppler emission profile or a 
complete redistribution in frequency) are inapplicable for describing the effect of neutral 
hydrogen on the Hell— >^HeI recombination kinetics and that the partial redistribution 
should be taken into account by using either numerical calculations (Switzer and Hirata 
2008), or fitting formulas (Wong et al. 2008), or an analytic solution (this paper), or a 

^Nh is the total number of hydrogen atoms and ions 
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combination of these approaches (Rubino-Martin et al. 2007). 

Our calculated values of N^/Nh agree with those calculated by Rubino-Martin et 
al. (2007) (for the model with a partial frequency redistribution of the resonance photons 
produced during transitions in the Hel singlet structure) with a relative accuracy of at least 
10"'^ (see Fig. [6l bottom panel). Our values of Ne/Nn calculated in the approximation of 
a complete redistribution of Hel resonance photons in frequency (i.e., using Eq. (l22l) for 
both 2^P l^S and 2^P l^S lines) agree with those of Rubino-Martin et al. (2007) 
(for a complete redistribution in frequency in the 2^P l^S and 2^P l^S lines) with 
a relative accuracy of at least 4 ■ 10~^ (see Fig. [6l bottom panel). 

Our additional results to be compared with those of other authors are the de- 
pendences of the probabilities Pfg = {Pfg + P}f) for the Hel 2^P l^S (Fig. W) and 
2^P ^1^5* (Fig. [HI) transitions on redshift z for various models of the redistribution of 
Hel resonance photons in frequency upon their scattering in lines. The calculated values 
of these quantities are in satisfactory agreement with those of Switzer and Hirata (2008) 
and Rubino-Martin et al. (2007). 

In conclusion, we calculated N^/Nh for various values of the cosmological param- 
eters. It should be noted that the z dependence of N^/Nh changes only slightly when 
varying the fraction of nonrelativistic matter VLm within the range 0.24 — 0.30 and the 
Hubble constant Hq within the range 65 — 75 km-s~^Mpc~^ (in practical calculations, for 
example, for the analysis of the CMB anisotropy spectrum, these changes may be ignored 
at modern level of experimental data; these results are not presented here graphically). 
The values of N^/Nh calculated by varying the fraction of baryonic matter Vt}, are pre- 
sented in Fig. [9] (top panel). N^/Nh is a monotonic function of fib- The Hell— >HeI 
recombination occurs at earlier epochs as fib increases. Figure [9] (bottom panel) presents 
the relative difference between the value of N^/Nh calculated using our model and the 
value of Ne/Nn calculated using the model by Wong et al.(2008) (i.e., in the Wong-Moss- 
Scott approximation). We see from [9] that the results of the calculations based on our 
model and the model by Wong et al. (2008) agree with a relative accuracy of at least 
2.2 • 10~^ (in the number density of free electrons) as the baryonic matter density varies 
in the range fib = 0.02 - 0.06. 

6 Conclusions 

We calculated the Hell— *^HeI recombination kinetics by including the effect of neutral 
hydrogen. We additionally took into account the partial redistribution of Hel resonance 
photons in frequency upon their scattering in the Hel 2^P —> l^S line (Eqs. 1351- 1441). 

It is shown that the calculated relative numbers of free electrons N^/Nh could be 
satisfactorily reconciled with the results of recent studies of the Hell— *^HeI recombination 
kinetics (Rubino-Martin et al. 2007; Wong et al. 2008) using the formulas derived here. 
The achieved accuracy is high enough for the Hell— >^HeI recombination model suggested 
here to be used in analyzing the CMB anisotropy data from future experiments (Planck 
and others). This accuracy is also high enough to calculate the intensities, frequencies, 
and profiles of the Hel recombination lines formed during the cosmological Hell— >^HeI 
recombination (Dubrovich and Stolyarov 1997; KhIV07; Rubino-Martin et al. 2007). 

It should also be noted that using Eqs. ((201), ((221), ((23]), ((351), ([39] - [44]) derived 



16 



here, along with Eq. ([29l) derived by Grachev (1988), we can find not only the dependence 
of on 7, but also the dependences of on the Voigt parameter a and the fraction 
of coherent scatterings A (and, accordingly, the fraction of incoherent scatterings (1 — A), 
if the atomic transitions produced by electron collisions may be neglected). This will be 
important in further detailed theoretical studies of the Hell— i>HeI recombination. 
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Figure 1: Model energy level diagram for the Hel atom used here. 
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Figure 2: Ratio of the Hel and HI absorption coefficients at the central frequency of the 
/ — > (? line (7) versus redshift z for the Hel 2^P — (solid curve) and Hel 2^P — 
(dashed curve) transitions. The dotted straight lines indicate the unit and 71 = 5 • 10^ 
levels. The value of 71 = 5 • 10^ is explained in the text. 



21 



4 • 



Oh 



10 = 



0,02 - 

0,01 K"' 

< 



-0,01 



1 — I I I 1 1 1 II 1 — I I I 1 1 III 1 — I I I 1 1 1 II 1 — I I I 1 1 1 II 1 — I I I 1 1 1 1 




I I M il l 



J I I I I I I I I 



10 



10^ 



10 10 
Y 



10' 



10 



10 



Figure 3: Top panel: Destruction probability of a Hel 2^P resonance photon as it 

interacts with neutral hydrogen versus 7 for various models of absorption and redistri- 
bution in frequency in the Hel 2^P line: the dash-dotted curve corresponds to the 
use of the Doppler absorption profile, the solid curve corresponds to a partial redistribution 
in frequency (Eq. (l35l) ). the dashed curve represents the Wong-Moss-Scott approxima- 
tion, and the dashed curve with two dots corresponds to the use of the Voigt absorption 
profile for a complete redistribution in frequency. The Voigt parameter is a = 1.7 ■ 10^^, 
and the single-scattering albedo is defined by the relation (1 — A) = 2.5 ■ 10~^. Bottom 
panel: The relative error in when calculated from the approximate formulas f l22l ) (the 
dashed curve with two dots corresponding to the use of the Voigt absorption profile for a 
complete redistribution in frequency) and (l20l) (the dash-dotted curve corresponding to 
the use of the Doppler absorption profile) compared to the numerical calculation based 
on the exact formula ffTSl). 
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Figure 4: Top panel: Destruction probability of a Hel 2^P l^S resonance photon as 
it interacts with neutral hydrogen versus 7 for various absorption profiles in the Hel 
2^P l^S line: the dashed curve corresponds to the use of the Doppler absorption profile 
and the solid curve corresponds to the use of the Voigt absorption profile for a complete 
redistribution in frequency. The Voigt parameter is a = 10~^ and the single-scattering 
albedo is approximately equal to zero A ~ (typical values for the Hel 2^P l^S 
transition at the Hell— *^HeI recombination epoch). Bottom panel: The relative error in 
P^ when calculated from the approximate formulas (1221) (the solid curve corresponding to 
the use of the Voigt absorption profile for a complete redistribution in frequency) and (l20l) 
(the dashed curve corresponding to the use of the Doppler absorption profile) compared 
to the numerical calculation based on the exact formula (fTSl). 
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Figure 5: Top panel: Destruction probability of a Hel 2^P l^S resonance photon 
as it interacts with neutral hydrogen Pq (l29ll versus 7 (solid curve), asymptotics (l33l) 
corresponding to 7 -C 71 (dashed curve) (Chugai 1987; see the text), and asymptotics 
([34II corresponding to 7 71 (dash-dotted curve). The Voigt parameters is a = 1.7- 10~^ 
and the single-scattering albedo is defined by the relation (1 — A) = 2.5-10"^ (typical values 
for the Hel 2^P l^S transition at the Hell— >HeI recombination epoch). Bottom panel: 
The relative error in Pq approximately calculated from Eqs. fl30l) and ([311) compared to 
the calculations based on Eq. (l29l) . 
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Figure 6: Top panel: Relative number of free electrons N^/Nh versus redshift z: the 
dashed curve represents the equilibrium recombination, the dotted curve represents the 
recombination according to the model without any absorption of Hel resonance photons 
by neutral hydrogen; the dash-dotted curve represents the recombination according to 
the model by KhlVOT; the solid curve represents the recombination according to this 
work (Eqs. ^ and ^ describe the Hel 2^P l^S and Hel 2^P l^S transition 
probabilities, respectively); the dashed curve with two dots represents the recombina- 
tion calculated using the approximation of a complete redistribution in frequency in Hel 
resonance lines (Eq. ( l22i ) describes the probability of both Hel 2^P — > l^S and Hel 
2^P l^S transitions). All curves were calculated for the cosmological parameters 
adopted in Rubino-Martin et al. (2007). Bottom panel: The relative difference in the 
fraction of free electrons for different recombination models: the solid curve indicates the 
difference between our results and the results by Rubino-Martin et al. (2007) for the mod- 
els with a partial redistribution in frequency; the dashed curve with two dots indicates 
the difference between our results and the results of Rubino-Martin et al. (2007) for the 
models with a complete redistribution in frequency; the dash-dash-dot curve indicates the 
difference between the results based on the model by Wong et al. (2008) and those by 
Rubino-Martin et al. (2007) for the model with a partial redistribution in frequency. 
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Figure 7: Probability of the uncompensated Hel 2^P — > l^S transitions Pbg versus redshift 
z for various models of absorption and redistribution in frequency in Hel 2^P — > l^S line: 
the dash-dotted curve corresponds to the use of the Doppler absorption profile, the solid 
curve represents a partial redistribution in frequency (this paper), the dash-dash-dot curve 
represents a partial redistribution in frequency (Rubino-Martin et al. 2007), the dashed 
curve represents the Wong-Moss-Scott approximation, the dashed curve with two dots 
corresponds to the use of the Voigt absorption profile for a complete redistribution in 
frequency, and the dotted curve represents the Sobolev probability Pj^ (calculated from 
Eq. (fT3l) ). The curves were calculated for the cosmological parameters adopted in Rubino- 
Martin et al. (2007). 
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Figure 8: Probability of the uncompensated Hel 2^P — > l^S" transitions Pb'g versus redshift 
z: the sohd curve represents our result (Eq. (l22l) was used for the calculations, which 
corresponds to the use of the Voigt absorption profile for a complete redistribution in 
frequency), the dash-dash-dot curve represents the result by Rubino-Martin et al. (2007), 
and the dotted curve represents the Sobolev probability P^g (calculated from Eq. f|T3ll). 
The curves were calculated for the cosmological parameters adopted in Rubino-Martin et 
al. (2007). 
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Figure 9: Top panel: Relative number of free electrons N^/Nh versus redshift z for 
various baryonic matter densities: the dashed, solid, and dash-dotted curves correspond 
to VLb = 0.02, Qb = 0.04, and fib = 0.06, respectively. The remaining parameters of the 
cosmological model are given in Table El Bottom panel: The relative difference in the 
number density of free electrons between the calculations using our model and those using 
the model by Wong et al. (2008) (the Wong-Moss-Scott approximation): the dashed, solid, 
and dash-dotted curves correspond to Qb = 0.02, Qb = 0.04, and Qb = 0.06, respectively. 
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